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OPTIMAL  RADAR  DOPPLER  PROCESSORS 


1.  INTRODUCTION 

Many  receivers  in  communications,  radar,  and  sonar  are  essentially  finite-memory 
digital  filters.  In  its  canonical  form,  a  finite-memory  linear  filter  is  equivalent  to  a  trans¬ 
versal  filter,  which  is  a  tapped  delay  line  whose  tap  outputs  are  weighted  and  summed. 
Therefore,  optimization  procedures  for  transversal  filters  would  have  wide  use  in  com¬ 
munications  and  detection  systems. 

One  such  application  of  the  transversal  filter  is  detection  of  moving  targets  by  radar 
systems  in  ground-clutter  or  sea-clutter  backgrounds  (Radar  returns  from  stationary  objects 
are  commonly  referred  to  as  clutter).  This  function  is  usually  performed  by  the  moving- 
target  indicator  (MTI).  The  MTI,  which  is  a  special  transversal  filter,  detects  moving  tar¬ 
gets  by  sampling  the  doppler  shifts  of  returns  from  fixed  objects  and  moving  targets  (1). 
lhe  moving  targets  produce  a  doppler  shift  and  are  passed  by  the  MTI;  returns  from  fixed 
objects  are  filtered  out.  Theoretically  the  doppler  filters  used  by  MTI  radars  are  relatively 
simple  in  that  they  process  a  small  number  of  returns,  typically  two  to  four.  If  modem 
digital  filtering  technology  (2)  is  used  many  more  returns  can  be  processed.  Optimization 
procedures  for  these  more  complex  filters  are  needed  to  account  for  the  effects  of  un¬ 
known  signal  and  clutter  parameters.  Because  of  the  random  nature  of  many  types  of 
clutter  returns,  in  this  research  clutter  is  treated  as  nonwhite  noise. 

Most  of  the  work  to  date  involving  moving-target  detection  has  been  concerned  with 
developing  a  theoretical  optimum  and  then  showing  that  the  performance  of  the  simple 
doppler  filters  in  use  is  very  close  to  the  optimum.  R.  C.  Emerson  (3)  has  developed  a 
method  for  minimizing  the  response  of  a  filter  to  clutter.  An  optimization  procedure 
developed  by  S.  P.  Applebaum  (4)  maximizes  the  signal-to-clutter  ratio  at  the  output  of 
such  a  filter  when  the  doppler  shift  of  the  signal  is  known. 

The  method  introduced  in  this  report  is  based  on  an  extension  of  the  Applebaum 
procedure  and  maximizes  the  signal-to-clutter  ratio  at  the  output  of  the  filter  when  the 
frequency  or  the  doppler  shift  of  the  signal  is  unknown  and  the  noise  is  nonwhite.  The 
investigation  starts  with  a  mathematical  development  of  the  optimization  equations,  using 
the  maximum-lif.elihood-ratio  test.  The  optimal  receiver  structure  depends  on  the 
covariance  functions  of  the  signal  and  the  noise,  which  are  not  usually  known  a  priori. 

In  the  second  part  of  the  investigation,  this  optimization  procedure  is  applied  to  the  MTI 
problem  to  illustrate  the  improvements  over  conventional  MTI.  A  generalized  doppler 
processor  is  developed  by  dividing  the  doppler  space  into  regions  and  optimizing  a 
processor  for  each  region.  The  results  of  these  processors  are  compared  with  those  of 
conventional  coherent  integration  filters  (5). 
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2.  OPTIMIZATION  CRITERION 

The  maximum-likelihood-ratio  receiver  processes  the  received  data  in  such  a  way  that 

ratio 

^  _  P(x!x  =  s  +  n) 

P(x|x  =  n) 

is  maximum.  Where  the  received  data  vector  is 
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P(x|x  =  *  +  n)  is  the  conditional  probability  of  receiving  vector  x  when  signal  s  is  trans¬ 
mitted.  P(xlx  =  n)  is  the  conditional  probability  of  receiving  vector  x  when  no  signal  is 
present.  When  these  probability  distributions  are  Gaussian  (that  is  when  noise  n  has  a 
Gaussian  distribution),  maximizing  the  likelihood  ratio  corresponds  to  maximizing  the 
signal-to-noise  ratio. 

The  optimal  transversal  filter  is  one  that  maximizes  the  output  signal-to-noise  ratio. 
The  output  signal  and  output  noise  refer  to  the  squares  at  the  absolute  values  of  the  signal 
and  the  noise  respectively.  The  expected  value  of  the  output  signal  is 

Ps  =  ajMga*,  (1) 


where 

Ms  is  the  signal  covariance  matrix 

a  is  the  weight  vector 

T  indicates  the  transpose  of  the  vector 

(*)  indicates  the  complex  conjugate  of  the  vector. 

The  expected  value  of  the  output  noise  is 

PN  =  aTAfva* 

where  MN  is  the  noise  covariance  matrix.  When  the  expected  value  of  the  input  signal 
and  the  expected  value  of  the  input  noise  are  normalized  to  unity,  the  improvement  in 
signed -to-noise  ratio  for  a  transversal  filter  can  be  defined  as 
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where 


jc.  =  exp  L l2nfd{i  -  1)7']. 


(3) 


Equation  (3)  is  the  ith  sample  of  the  input  function  of  time  defined  by  Eq.  (A17).  It  is 
complex  because  Eq.  (A17)  is  complex  i.e.,  it  is  a  function  not  of  a  real  frequency  but  of 
a  frequency  shift.  With  the  demodulation  procedure  outlined  in  Appendix  A,  negative  as 
well  as  positive  frequency  shifts  can  be  detected.  Therefore,  if  the  carrier  frequency  is 
removed  and  the  doppler  shift  is  considered  as  a  real  frequency,  complex  signals  and  com¬ 
plex  autocorrelation  functions  will  result  in  practical  applications.  A  further  result  of  this 
is  that  if  the  mean  square  of  the  absolute  value  of  this  signal  is  considered  as  the  power 
spectrum  this  "power  spectrum"  may  not  he  symmetrical  with  respect  to  the  "zero- 
frequency”  axis. 

Assuming  stationarity  and  ergodicity,  we  find  the  element  (i.  k)  of  the  covariance 
matrix  by  taking  the  expected  value  with  respect  to  the  unknown  doppler  shift. 


E 


xixk 


exp  [2jirfd{i  -  k)T]  P[fd )dfd 


where  Pifd )  is  given  above. 

It  follows  that 

....  ,  ,,  sin  ir(i  -  k) 

mt*  3  «P  n{i  .  <4> 

for  i  =  1,  2,  N  and  k  =  1,  2,  ....  IV. 

From  Eq.  (4), 

1,  i  =  k 

mi.k  =  ‘ 

C,  i*k. 

Therefore,  Ihe  signal  covariance  matrix  for  unknown  doppler  shift  is 

Ms  =  /  (5) 
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where  /  is  the  identity  matrix, 
into  Eq.  (2): 


TTie  improvement  factor  is  found  by  substituting  Eq.  (5) 


ara 


* tin 


*TMN  a 


(6) 


Maximizing  Eq.  (6)  is  equivalent  to  minimizing 

Py  a  3  yvAfy  a  t 

where  Pv  is  the  output  noise  power.  This  has  been  shown  (3)  to  be  minimum  when 
weight  vector  a  is  chosen  to  be  the  eigenvector  that  results  in  the  smallest  eigenvalue  of 
My,,.  The  improvement  in  signal-to-noise  ratio  is  given  by  the  reciprocal  of  this  eigenvalue. 


2.2.  Multiple  Filters 

It  has  been  shown  (5|  that  N  independent  filters  can  be  generated  with  N  samples  of 
the  input  data.  By  dividing  the  rr.alyzing  bandwidth  (0.  1  IT)  into  N  equal  intervals,  a 
filter  can  be  designed  to  detect  optimally  a  signal  in  each  interval.  In  this  way,  the  entire 
analyzing  bandwidth  can  be  covered,  but  the  signal  for  a  particular  filter  can  be  assumed 
to  have  equal  probability  of  occurring  anywhere  within  the  interval  covered  by  that  filter 
instead  of  anywhere  within  the  analyzing  bandwidth.  This  additional  information  on  the 
signal  can  be  used  to  improve  the  detectability  of  the  signal.  The  goal  is  to  find  the 
weight  vector  to  maximize  Eq.  (2)  for  a  signal  whose  doppler  shift  is  known  to  fall  into 
an  interval  ((2m  -  l)/2/V7\  (2m  +  DI2NT),  m  =  0,  1,  . . . ,  N  -  1. 


Derivation  of  Optimal  Weight  Vector— The  optimal  weights  are  given  by  the  vector 
A  that  maximizes  Eq.  (2): 


*$fn 


aT^Sa 

aTAf*a* 


To  find  the  desired  vector,  the  following  definitions  will  be  used: 


M's  -  WtMsW* 

(7) 

M'n  - 

(8) 

Further,  matrix  W  will  be  defined  in  such  a  way  that 

(9) 

Rewriting  Eq.  (7)  yields 

Ms  “(R^r’A/jtlF*)-1. 

(10) 

Rewriting  Eq.  (8)  and  making  use  of  Eq.  (9)  results  is 

Mv  =  (lVr)-,(W'*r1. 

(11) 
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Substituting  Ews.  (10)  and  (11)  into  Eq.  (2),  we  have 


aT(WT)-,A/;.(W’)  'a* 

**ln  •T{wTr'iw'r'f 


Define  a  vector  f  as 


Then 


which  can  be  written  as 


f  =  tr*a. 


sin 


W* 

fTt* 


-  w 

llfll2-frf*  =  1. 


(12) 


(13) 

(14) 


Equation  (13)  is  maximized  when  f  is  chosen  to  be  the  eigenvector  that  results  in  the 
largest  eigenvalue  of  M's. 


Thus,  it  has  been  shown  that  the  optimum  weight,  when  the  doppier  shift  is  known 
to  be  within  some  interval  of  values,  is 


aopt  =  **  <15> 

where  w  is  defined  by 

w*w  (16) 

and  f  is  the  eigenvector  that  produces  the  largest  eigenvalue  of  w^AfjW*. 


Signal  Covariance  Matrix— The  signal  covariance  matrix  can  be  derived  in  a  way 
similar  to  that  used  when  the  doppier  shift  is  assumed  to  be  completely  unknown.  The 
difference  is  that  for  this  case  the  signal  is  assumed  to  have  equal  probability  of  occurring 
anywhere  within  a  region  ((2 n  -  l)/2 NT,  (2 n  +  1)I2NT),  n  =  0,1,  N  -  1,  where  T  is 
the  delay  of  each  delay  line  and  n  is  an  index  used  to  select  the  region  of  the  analyzing 
bandwidth  (0,  1/T)  for  this  filter.  Therefore,  the  probability  density  function  for  the 
doppier  frequency  fd  is 


2n-l  2/1  +  1 

2 NT  2 NT 

otherwise. 
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This  leads  to  the  element  (i,  k )  of  the  signal  covariance  matrix. 


m 


.>  =  «P 


l2irn(i  -  k) 


.1  sin 


V 


n(i-k) 


NT 


Mi  -k) 


NT 


(17) 


fori  =  1,2 . N 

fe  =  1,  2 . N 

n  =  0,  1 . N  -  1. 

If  Eq.  (17)  is  used,  the  signal  covariance  matrix  for  a  signal  whose  doppler  shift  is 
known  to  lie  within  the  interval  ((2 r  -  l)/2 NT,  (2/i  +  1)/2jV7")  can  be  generated.  Then, 
Eqs.  (15)  and  (16),  we  can  compute  the  optimal  weight,  vector.  For  the  case  where  A/:1 
does  not  exist,  the  determinant  of  .Wv  is  zero,  which  implies  that  at  least  one  eigenvalue 
of  A/v  is  zero.  The  optimum  weight  vector  for  this  case  is  the  eigenvector  that  produces 
one  of  the  zero  eigenvalues  of  M . 


3.  OPTIMAL  RADAR  MTI  PROCESSORS 

The  optimization  procedures  developed  in  the  preceding  section  can  be  applied  to 
the  detection  of  moving  targets  by  a  radar  system.  The  doppler  shift  of  the  returns  from 
moving  targets  are  unknown  in  general.  Therefore,  they  are  usually  assumed  to  have  equal 
probability  of  occurring  anywhere.  For  this  reason,  the  detection  system  must  be 
optimized  to  detect  a  target  that  has  any  doppler  shift  within  the  analyzing  bandwidth  of 
the  radar  (0,  1/T),  where  1  /T  is  the  pulse  repetition  frequency  iPRF). 


3.1.  Conventional  MTI 

The  conventional  MTI  canceler  delays  the  returns  of  a  given  transmittev  ,  uise  and 
subtracts  them  from  the  returns  of  the  next  transmitted  pulse.  A  number  of  cancelers  n 
cascaded  as  shown  in  Fig.  2  are  equivalent  to  an  (n  +  l)-sample  transversal  filter  with 
weights  corresponding  to  the  nth-degree  binomial  coefficients  with  alternating  signs.  The 
power  transfer  functions  for  these  filters  can  be  computed  using  P  =  aA/a  with  weight 
vector  a  given  by  the  binomial  coefficients  with  alternating  signs  and  the  covariance  matrix 
given  by 

mj  k  =  cos  2nf0(k  -  i)T. 

7.  These  curves  are  normalized  by 


The  results  are  shown  in  Fig.  3  for  n  =  1  through  n 
dividing  the  gain  by  the  maximum  gain 


n*  1 


L 


(18) 


where  a(  are  the  binomial  weights. 
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Ki||.  2  —  An  n-slage,  or  n-canceler,  MTI 


Fig.  3  — Normalized  power  transfer  function  P,  for  an  MTI  having  the 
indicated  number  of  cancelers 


The  improvement  in  signal-to-noise  ratio  can  be  computed  using  Eq.  (2).  Weight 
vector  a  represents  the  binomial  weights.  Sine®  the  doppler  shift  is  considered  to  have 
equal  probability  of  assuming  any  value,  signal  covariance  matrix  Afs.  was  shown  to  be  the 
indentity  matrix.  It  remains  only  to  define  the  clutter  (noise)  covariance  matrix  Afv  in 
order  to  evaluate  Eq.  (2). 

In  a  typical  clutter  environment,  the  clutter  energy  received  by  the  radar  is  made  up 
of  the  returns  from  a  large  number  of  scatterers  within  a  resolution  cell  of  the  radar. 

This  resolution  cell  is  determined  by  the  antenna  pattern  and  the  transmitted  pulse  width. 
The  scatterers  are  randomly  distributed  within  the  resolution  cell,  and  they  generally  have 
random  internal  motion  such  as  the  fluttering  of  leaves  on  trees  or  falling  raindrops. 
Therefore,  the  clutter  must  be  described  as  a  random  variable. 

To  design  an  optimum  detection  system,  both  the  probability  distribution  of  the  clut¬ 
ter  amplitude  statistics  and  the  clutter  spectral  shape  must  be  known.  Nathanson  and  Reilly 
(6)  have  considered  the  effects  of  clutter  statistics  on  radar  performance  and  have  shown 
that  current  knowledge  does  not  allow  an  optimal  detection  criterion  to  be  specified  a  priori. 
With  this  dilemma,  one  is  naturally  led  to  adaptive  techniques  in  which  critical  clutter  param¬ 
eters  are  estimated  and  the  receiver  characteristics  are  adjusted  accordingly. 
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When  the  clutter  return  in  each  resolution  cell  is  due  to  many  statistically  independ¬ 
ent  scatterers,  it  follows  from  the  central  limit  theorem  that  this  clutter  can  be  described 
by  a  Gaussian  amplitude  distribution.  This  model  is  representative  of  many  types  of 
clutter,  such  as  rain,  forests,  and  sea  returns,  and  is  assumed  to  apply  for  this  optimization 
The  distribution  of  the  doppler  shift  of  clutter  has  been  found  to  be  highly  dependent  on 
the  type  of  clutter  and  on  weather  conditions,  particularly  wind. 

In  prior  analyses  of  MT1  systems,  a  Gaussian  distribution  of  the  clutter  doppler  shift 
was  used  to  evaluate  the  performance  of  the  MTI  system.  For  this  reason,  a  Gaussian 
distribution  will  be  used  to  compare  the  performance  of  MTI  systems  using  prior  theory 
with  the  performance  of  the  optimal  processor  developed  as  a  part  of  this  research.  The 
model  for  this  doppler  shift  is 


Pr<4>'C0exP 


-<4  ~Pf> 

2  ol 


21 


(19) 


where  u  and  of  represent  the  mean  and  variance  of  this  distribution  and  C0  is  the  energy 
level.  Tne  Foruier  transform  of  Eq.  (19)  is 

v(r)  =  ocCoV^2iexp(-2JT2o;r2 ->2^7).  (20) 


Since  the  receiver  design  is  not  affected  by  CQ,  let 

1 


C0  = 


Therefore,  element  ( i.  k)  of  ,W  is 

mi,k  =  e*Pl-  -  k)2r2  ->2w#Jtc(t  -  k)T] .  (21) 

With  Eq.  (21),  all  the  terms  of  Eq.  (2)  are  defined,  and  the  improvement  factor  for  a 
conventional  MTI  with  binomial  weights  can  be  computed.  The  results  are  shown  in  Fig. 

4  for  =  0  and  0  001  <  ofT  <  0.1  where  1  IT  is  the  PRF.  The  mean  has  been  set  equal 
to  zero  for  this  comparison,  although  many  types  of  clutter  do  not  have  a  zero-mean 
doppler  shift,  in  particular  clutter  from  clouds  and  rainfall  or  that  received  when  the  radar 
is  on  a  moving  platform.  If  this  nonzero  mean  is  not  taken  into  account,  a  smaller  im¬ 
provement  factor  results.  These  effects  will  be  considered  further  in  a  later  report. 


3.2.  Optimal  MTI 

The  procedure  for  deriving  the  optimal  MTI  weights  was  developed  in  Sec.  2.  It  was 
shown  that  when  any  value  for  the  doppler  shift  is  equally  probably,  the  optimum  weights 
are  given  by  the  eigenvector  that  produces  the  smallest  eigenvalue  of  MjS.  The  improve¬ 
ment  factor  is  the  reciprocal  of  that  eigenvalue. 

If  Eq.  (21)  is  used  to  generate  MN,  the  improvement  factor  for  the  optimal  MTI  is 
shown  in  Fig.  5  for  N  =  1  through  N  *  4,  where  N  is  the  number  of  cancelers.  For 
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comparison,  the  dotted  curves  represent  the  improvement  factor  for  the  conventional 
binomial  MT1.  The  corresponding  eigenvectors  (optimum  weights)  are  tabulated  in  Ap¬ 
pendix  B  for  iV  =  2  through  4.  A  single  canceler  (N  =  1)  is  optimized  with  weights  (1,  -1) 
for  all  values  of  ot.T.  These  are  the  binomial  coefficients  with  alternating  signs  for  N  =  1. 
In  general,  the  binomial  weights  are 

*,  =  <-! )'<*),  (  =  0,1,  ...,1V.  (22) 

For  a  double  canceler  (N  =  2),  the  binomial  weights  ran  be  normalized  to 

a  =  (-0.5, 1.0,  -0.5)=  i  (-  1,  2,  -  1). 

The  normalized  binomial  weights  for  N  -  3  and  4  are 

a  =  (-  .333, 1.0,  -  1.0,  0.333)  =  ~  (-  1,  3,  -  3, 1) 

3 

and 

a  =  <0.1667,  -  0.6667, 1.0,  -  0.6667,  0.1667)  -  ~  (1,  -  4,  6,  -4,1). 

6 

If  the  clutter  spectrum  is  very  narrow  (i.e.,  ofT  -*  0),  the  optimum  weights  are  very 
nearly  equal  to  the  binomial  weights.*  However,  the  small  differences  between  these 
optimum  weights  and  the  binomial  weights  result  in  an  appreciable  increase  of  the  im¬ 
provement  factor.  From  this  it  can  be  concluded  that  the  accuracy  of  the  weights  be¬ 
comes  very  critical  if  the  achievable  improvement  factor  is  large. 

The  transition  of  the  optimum  weights  as  the  spectral  width  of  the  clutter  doppler 
shift  increases  is  shown  in  Fig.  6  for  N  =  3.  This  is  typical  of  the  results  for  other  values 
of  N.  The  absolute  value  of  the  weights  toward  the  end  of  the  transversal  filter  (flj  and 
a.-),  in  general,  increases  as  the  width  of  the  clutter  spectrum  increases.  This  results  in  a 
filter  transfer  characteristic  that  has  a  narrower  main  lobe  and  higher  side  lobes,  as  shown 
in  Fig.  7.  The  dotted  curve  is  the  binomial  weighted  filter  with  no  sidelobes,  which  is 
shown  for  comparison  with  the  optimal  filter  for  a  wide  clutter  spectrum  (°CT  -  0.5). 

The  appearance  of  side  lobes  for  these  filters  is  revealed  by  careful  examination  of  the 
filter  characteristics  near  zero  doppler  shift. 

When  the  signal  doppler  shift  is  considered  to  have  equal  probability  of  any  value,  it 
has  been  shown  that  the  signal  covariance  matrix  is  the  identity  matrix.  As  the  clutter 
spectrum  width  ac  is  increased,  the  clutter  approaches  the  characteristics  of  white  noise 
in  that  the  clutter  covariance  matrix  approaches  the  identify  matrix.  For  this  case,  it  is 
seen  from  Eq.  (2)  that  the  improvement  factor  approaches  unity  (0  dB)  regardless  of  the 
weights.  Therefore,  the  accuracy  of  the  optimum  weight  is  less  critical  for  a  very  wide 
clutter  spectrum  and  a  small  maximum  achievable  improvement  factor.  This  is  shown  in 
Fig.  8,  where  the  clutter  spectral  width  of  T  is  increased  to  0.5. 

In  summary,  this  optimization  procedure  results  in  a  significant  additional  improve¬ 
ment  in  signal-to-clutter  ratio  when  (a)  more  than  two  pulses  (N  >  1)  are  processed,  (b) 
the  clutter  spectrum  is  narrow,  and  (c)  the  filter  weights  are  accurate,  as  tabulated  in 


•This  can  he  seen  in  Appendix  B,  tables  Bl,  B2,  and  B3,  for  small  values  of  ncT. 
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Fig.  8  —  MTI  improvement  factor,  optimum  weight*  solid  curves, 
binomial  weights  broken  curves.  N  >  number  of  cancelert. 


Appendix  B.  If  the  standard  deviation  of  the  clutter  spectrum  is  0.01  times  the  PRF, 

Fig.  5  shows  that  the  additional  improvement  factor  achievable  with  optimum  weights  is 
about  2  dB  for  a  double  canceler  (N  =  2),  3.5  dB  for  N  =  3,  and  5  dB  for  N  *  4.  Al¬ 
though  the  weights  must  be  accurate,  Fig.  6  shows  that  they  change  very  slowly  as  the 
width  of  the  clutter  spectrum  is  changed.  Therefore,  if  the  approximate  width  of  the 
clutter  spectrum  is  known  for  a  given  application,  the  optimum  weights  can  be  computed 
and  adaptation  can  be  avoided. 


4.  GENERALIZED  JV-PORT  DOPPLER  PROCESSORS 

In  the  preceding  section  an  MTI  was  optimized  for  detecting  a  signal  whose  doppler 
shift  is  given  equal  probability  of  having  any  value.  It  was  shown  (Fig.  8)  that  using  the 
“best”  weights  gives  significantly  better  results  than  using  binomial  weights  if  the  returns 
from  a  large  number  of  pulses  are  processed.  However,  when  the  clutter  spectrum  be¬ 
comes  very  wide,  the  improvement  in  sign al-to -clutter  ratio  of  the  optimal  MTI,  as  well 
as  the  binomial  MTI,  approaches  0  dB,  as  shown  in  Fig.  8. 

The  improvement  in  signal-to-noise  ratio  of  both  the  optimal  and  the  binomial  MTI 
filters  with  a  signal-plus- white-noise  input  is  also  0  dB.  This  can  be  seen  by  referring  to 
Eq.  (2)  and  remembering  that  the  covariance  matrix  of  white  noise  is  the  identity  matrix, 
and  the  covariance  matrix  of  a  signal  whose  doppler  shift  is  unknown  is  also  the  identity 
matrix.  Therefore,  the  improvement  factor  is  unity  (0  dB)  regardless  of  the  weights  used. 

In  a  realistic  radar  system,  the  receiver  must  contend  with  both  clutter  and  white 
noise  produced  in  the  input  stage  of  the  receiver.  To  cope  with  clutter-plus-noise  inter¬ 
ference,  the  doppler  filters  are  generalized  to  form  N  independent  doppler  filters  from  N 
samples  of  the  input  data.  Each  filter  covers  a  fraction  (1  /NT)  of  the  analyzing  bandwidth 
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(0,  1  IT).  In  this  way  the  entire  spectrum  is  covered,  but  each  filter  is  optimized  only  for 
the  region  to  which  it  is  assigned. 


4.1.  Interference  Model  for  Clutter  Plus  Noise 

To  derive  a  model  for  clutter  plus  noise,  it  is  assumed  that  (a)  the  total  energy  of  the 
interference,  which  includes  clutter  plus  noise,  is  normalized  to  unity,  (b)  the  clutter  has  a 
Gaussian  spectrum,  (c)  the  noise  has  a  “white"  frequency  spectrum,  and  (d)  the  clutter 
and  noise  both  have  Gaussian  amplitude  probability  density  functions. 

Let  Ec  represent  the  power  of  the  Gaussian  spectrum  and  Ew  represent  the  power 
of  the  “white”  spectrum.  Then 


Eu  +  Ew  -  1 


and  the  total  interference  spectral  density  is 


P,(f) 


1 


w 


OcV^ 


exp< 


2  at 


>*EWT 


ufn-uff-jr) 


where  and  ac  are  the  mean  and  variance  of  the  Gaussian  spectrum  and 

fo,/<  0 

U(f)  M 

U,/>o. 

TTie  total  energy  of  the  interference  is  found  by  integrating  Eq.  (23): 


h 


(f)df~  1. 


The  Fourier  transform  of  Eq.  (23)  is 


4/(t)  »  (1  -  Ew )  exp(-  2n2o2r2  -;2ffpeT)  + 


w 


A, 

(?) 

The  element  (i,  k)  of  the  covariance  matrix  is 

mi  k  =  (1  -  Ew)  exp{-  2jt2ct2(j  -  k)2T2  -/2irjic(i  -  k)T] 

_  sin  rr (i  -  k)  ,vl 

+  EW  n{i„k)  °pW'-*)]. 


(?) 


(23) 


(24) 
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which  can  be  rewritten  as 

mik  -  (1  -Ew)exp[-  2v2o*(i -k)2T2  -/2jrpe(» -*)T]  +  Ew6(i-k)  (25) 

where 

1,  i-k 

Hi -it) 

[o,  i  *  k. 


4.2.  Effect  of  White  Noise  on  Optimal  MTI 

If  Eq.  (25)  is  used,  the  interference  covariance  matrix  can  be  generated.  If  Eq. 
(25)  is  compared  with  Eq.  (21),  it  can  be  seen  that  is 

M,  =■  (1  -  EW\MC  +  Ewl  (26) 

where  Mp  is  the  covariance  matrix  of  the  Gaussian  clutter  spectrum.  Identity  matrix  /  is 
the  covariance  matrix  of  the  white  spectrum.  It  has  been  shown  that  the  optimum  weight 
vector  for  an  MTI  is  the  eigenvector  that  gives  the  smallest  eigenvalue  of  M{.  To  find  this 
weight  vector,  the  eigenvalues  and  eigenvectors  of  Mt  will  be  calculated  from  the  eigen¬ 
values  and  eigenvectors  of  Mr. 

The  eigenvalue  equation  for  M}  is 


M,u  *  X,u 

where  u  is  an  eigenvector  of  Afj  and  is  an  eigenvalue.  If  Eq.  (26)  is  used, 

[(1  ~EW)MC  +£h,/]u*  X;u, 


which  can  be  rewritten  as 


where 


Mcu 


\  ~Ew 


1-E, 


w 


u-Xcu 


\ 


X,  -  Ew 


1-E, 


w 


(27) 


(28) 


(29) 


Since  Eq.  (28)  is  the  eigenvalue  equation  for  Afc,  it  follows  that  Mc  and  M{  have  the  same 
eigenvectors.  The  eigenvalues  are  related  by  Eq.  (29). 

The  eigenvalues  of  are  found  by  rewriting  Eq.  (29)  as 

X,  *  (1  -Ew)\.  +EW.  (30) 

When  Xc  is  a  minimum,  X,  is  also  a  minimum.  Therefore,  the  eigenvector  that  produces 
the  minimum  eigenvalue  of  Afc  also  produces  the  minimum  eigenvalue  of  Afj. 
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In  summary,  the  optimum  weight  vector  for  a  Gaussian  clutter  spectrum  is  also  the 
optimum  weight  vector  when  a  white  spectrum  is  added  to  the  Gaussian  spectrum.  The 
improvement  in  signal-to-noise  ratio  for  the  optimal  MTI  is  given  by  the  reciprocal  of  the 
minimum  eigenvalue  of  the  interference  covariance  matrix, 


& 


1 

<X/W 


where  l  '%  n  is  the  improvement 
Eq.  (30)  is  used. 


factor  for  a  Gaussian  spectrum  plus  a  white  spectrum. 


1  <!-*«,) 

77  ’  I  *Ew' 


‘tin 


tin 


If 

(31) 


where  lg/n  is  the  improvement  factor  for  a  Gaussian  spectrum  only.  Equation  (31)  shows 
that  the  improvement  factor  for  the  optimal  MTI  when  white  noise  is  added  is  related 
only  to  the  amount  of  white  noise  Ew  and  the  improvement  factor  for  a  Gaussian  spectrum 
1%/n-  eigenvectors,  or  optimum  weights,  are  not  affected.  This  means  that  the  deg¬ 
radation  of  both  the  conventional  and  optimal  MTI  by  white  noise  can  be  expressed  by 
Eq.  (31),  which  is  plotted  in  Fig.  9. 


Fig.  9  —  Effect  of  white  noise  on  the  improvement  factor  of 
both  the  optimum  and  the  binomial  MTI.  Is/n  ia  the  improve¬ 
ment  factor  if  the  interference  ia  clutter  alone.  rf/n  is  the 
improvement  factor  for  clutter  plus  noise.  is  the  fraction 
of  the  interference  that  ia  white  noise. 
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Figure  9  shows  that  when  lgjn  becomes  very  large,  /'/n  is  determined  (and  limited)  by 
Ew.  Under  this  condition,  the  achievable  improvement  in  signal-to-noise  ratio  cannot  be 
increased  by  processing  more  samples  (increasing  JV),  which  would  increase  Igjn  but  would 
not  increase  l'sjn . 

1.  is  also  increased  by  decreasing  clutter  spectrum  width  a,.  However,  when  Igjn 
becomes  very  large,  l'g,n  is  limited  by  the  white  noise  energy  k’w.  Hiis  is  shown  in  Fig. 

10  for  JV  =  2  through  JV  =  4. 


4.3.  Optimization  of  the  iV-Port  Doppler  Processor 

To  improve  the  performance  for  wide  clutter  spectra  and  white  noise,  a  contiguous 
bank  of  JV  filters  is  formed  by  applying  jV  independent  weight  vectors  to  the  JV  sample 
input.  Each  filter  is  optimized  to  cover  a  portion  of  the  analyzing  bandwidth.  The 
optimum  weights  are  specified  in  this  section  and  the  performance  of  these  filters  are 
illustrated. 

Hie  optimum  weight  vectors  for  signals  whose  doppler  shift  is  known  to  be*  within  a 
bounded  region  was  considered  in  the  second  section,  where  the  optimum  weights  are 

aopt  *  Wi-  (32) 

Hie  Matrix  W  is  defined  by 

WW*T  -  MjJ .  (33) 

Vector  f  is  the  eigenvector  associated  with  the  largest  eigenvalue  of  Af', 
where 


M's  =  W’tMsW.  (34) 

Performance  Against  Clutter  With  Gaussian  Spectrum— As  shown  above,  the  optimum 
weights  and  the  improvement  factor  can  be  found  if  the  covariance  matrix  of  clutter  A/, 
and  the  covariance  matrix  of  the  signal  M  are  known.  Equation  (21)  can  be  used  to 
generate  the  covariance  matrix  for  a  Gaussian  clutter  spectrum.  Equation  (17)  can  be 
used  to  generate  the  covariance  matrix  for  a  signal  whose  doppler  has  a  uniform  probability 
distribution  within  the  region  ((2 n  -  1)/2jV7',  (2 n  +  1)/2JV7),  n  =  0,  1,  ....  JV. 

If  Eqs.  (17)  and  (21)  are  used  to  generate  ,Vf  and  Mc ,  the  optimum  weights  and  the 
improvement  factor  are  found  for  JV  =  2  through  N  =  5.  The  results  are  shown  in  Fig.  11 
for  pc  *  0  and  0.005  <  a CT  <  0.5.  The  S  filters  are  shown  by  the  solid  curves  graph, 
the  filters  are  numbered  starting  with  the  one  at  zero  frequency.  The  average  improve¬ 
ment  factor  for  all  JV  filters  are  shown  by  the  dotted  curves.  Since  the  solid  curves 
represent  the  average  improvement  factor  for  a  doppler  shift  with  uniform  probability 
distribution  within  the  region  of  the  filter,  and  since  the  entire  analyzing  bandwidth  (0, 

1/T)  is  covered  by  the  JV  filters,  the  dotted  curves  represent  the  average  improvement 
factor  for  a  doppler  shift  with  a  uniform  distribution  over  the  analyzing  bandwidth. 
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(a)  Two-port  processor 


!  t  'M .  >  R!  iM 


:-'0r  ,pH  !  'H; 


(h)  Three-port  processor 

Fig.  1  1  -  Improvement  factor  of  an  /V-porl.  floppier  processor  compared  with 
that  of  an  optimum  /V-pulse  MTI. 
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(c)  Four-port  processor 


*  i  •  ■]  i  h  ''JTT  i  *.j'r .  Vc*  '  VifS  :'Ht 

(d)  Five-port  processor 

Fin.  11  (Continued)  —  Improvement  factor  of  an  N-port  doppler  processor 
compared  with  that  of  an  optimum  A/'-pulse  MTI. 
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For  comparison,  the  improvement  factors  using  the  optimum  weights  derived  for  the 
MTI  are  shown  by  dotted  curves  in  Fig.  11.  The  improvement  factor  for  the  optimal  MTI 
and  the  average  improvement  factor  for  the  Af-port  doppier  processor  coincide  until  the 
clutter  spectrum  becomes  very  wide.  The  reason  for  this  can  be  understood  by  examining 
Fig.  12,  which  shows  the  optimum  weights  of  each  niter  for  a  two-pulse  processor  ( N  =  2) 
The  □  curves  are  the  real  parts  of  the  weights  and  the  X  curves  are  the  innginary  parts. 
Both  filters  have  weights  of  (1,  -  1)  until  the  clutter  spectrum  becomes  ve"y  wide  («.  T  > 
0.15(.  Therefore,  since  both  filters  have  the  same  weights,  only  one  is  needed  for  a  T  < 
0.15.  That  one  filter  has  the  same  weights  as  the  optimum  MTI  for  N  -  2.  It  follows 
that  the  improvement  factor  would  be  the  same  in  eacn  case. 

When  acT  >  0.15,  the  weights  on  filter  1  become  (1,1).  At  this  point  the  two  fil¬ 
ters  are  different,  so  that  two  filters  are  needed  ard  the  improvement  factor  is  greater 
than  for  the  optimum  MTI,  as  was  seen  in  Fig.  11.  The  transfer  characteristics  of  the  two 
filters  are  shown  in  Fig.  13.  The  □  curve  represents  the  optimum  2-pulse  MTI.  The  A 
and  X  curves  represent  the  2-port  processor  for  acT  >  0.15. 

The  optimum  weights  for  the  four  filters  of  a  4-port  processor  are  shown  in  Figs. 
14-17.  When  the  clutter  spectrum  is  narrow  ( afT  <  0.05),  the  optimum  weights  of  all 
four  filters  are  identical  and  they  are  the  same  as  the  optimal  weights  of  the  optimal  4- 
pulse  MTI.  Therefore,  again,  only  one  filter  is  needed  for  narrow  clutter  spectra,  and 
that  filter  is  the  optimal  MTI. 

An  IV-point  discrete  Fourier  transform  (DFT)  is  given  by 


akn  =  «P 


(n  -  1)<*  -  1) 
N 


”,  k  1,  2 . iV. 


This  operation  corresponds  to  N  contiguous  filters.  For  N  =  2,  the  DFT  weights  for  the 
two  filters  f  t  and  f2  are 


f,  -11,1) 

and 

f2  -11,-1). 

These  weights  correspond  to  those  of  the  optimal  2-port  processor  for  ?  wide  clutter 
spectrum.  For  N  =  4,  the  DFT  weights  fer  the  four  filters  are 

f, -a,  i,i,D 

f2  -  o,  i.  -/,  -1) 

#3  -  (-1.1, -1.1) 

U  l,/,-l). 

These  weights  are  similar  to  those  of  the  optimal  4-port  processor  for  a  wide  clutter 
spectrum.  The  optimal  4-port  weights  are  modified  by  a  “window”  (0.83,  1,  1,  0.83). 

!t  is  well  known  that  the  optimal  window  is  uniform  (1,  1,  1,  1)  when  the  noise  is  white. 
When  the  clutter  spectrum  is  very  wide,  the  optimal  N-port  processor  is  very  close  to  an 
/V-point  DFT. 
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In  summary,  the  optimal  N -port  processor  is  a  single  filter  (i.e.,  the  optimal  MTI) 
when  the  clutter  spectrum  is  narrow  and  approaches  an  N-point  DFT  when  it  is  wide.  In 
the  transition  region,  the  improvement  factor  is  very  low  (Fig.  11),  and  the  weights  may 
change  drastically  for  small  changes  in  the  width  of  the  clutter  (Figs.  14-17).  However, 
the  filter  gain  changes  very  slowly  in  this  region  (Fig.  11),  and  this  implies  that  the  values 
of  the  weights  are  not  critical  in  this  region. 

The  filter  transfer  characteristics  of  three  of  the  four  filters  of  the  4-port  processor 
are  shown  in  Fig.  18.  The  □  curves  represent  the  filters  when  the  clutter  is  narrow 
( acT  *  0.005).  All  four  filters  are  the  same  for  this  condition  and  all  four  are  also  the 
same  as  the  4-pulse  optimal  MTI  (Fig.  3).  Not  all  of  the  filters  have  maxima  in  their 
detection  regions  when  the  clutter  is  narrow,  as  one  might  expect. 

The  X  curves  represent  the  optimal  filters  when  the  clutter  is  wide  ( acT  >  0.5). 

Notice  t!iat  each  is  very  close  to  a  (sin  x)x  shape,  the  shape  a  DFT  would  have.  The 
peak  of  each  filter  is  centered  in  the  detection  region  for  this  case. 

The  A  curves  represent  the  filter  shapes  in  transition  region  (acT  *  0.1).  The  peaks 
are  closer  to  the  detection  region  than  for  narrow  spectrum  clutter,  and  the  mainlobes  of 
the  filters  are  narrower  and  the  sidelobes  are  higher  than  the  (sin  x)x  shape.  This  cor¬ 
responds  to  an  inverse  taper  on  the  window  function;  i.e.,  the  end  samples  are  weighted 
heavier  than  the  center  samples. 

Similar  curves  of  the  filter  shapes  for  N  *  3  and  N  *  5  are  shown  in  Appendic  C. 

Performance  Again.it  Clutter  Plus  White  Noise—  It  was  shown  above  that  the  optimal 
MTI  does  not  improve  the  signal-to-white-noise  ratio.  Thus,  the  MTI  improvement  factor 
becomes  limited  by  the  white  noise  (Fig.  10),  and  processing  more  pulses  or  reducing  the 
clutter  spectral  width  ocT  does  not  lead  to  a  higher  improvement  factor.  This  difficulty 
is  overcome  by  the  <V-port  doppler  processor. 

If  Eq.  (25)  is  used  to  generate  interference  covariance  matrix  \J,  and  Eq.  (17)  to 
generate  signal  covariance  matrix  Ms,  one  can  find  the  largest  eigenvalues  of  M's  which, 
along  with  the  associated  eigenvectors,  give  the  improvement  factor  and  the  optimum 
weights  for  the  jV-port  doppler  processor.  The  results  of  this  computation  are  shown  in 
Figs.  19-22  for  2-  through  5-port  doppler  processors.  The  improvement  factors  of  the  N 
ports  are  averaged  and  shown  along  with  the  improvement  factor  of  the  optima]  N-pulse 
MTI  for  comparison. 

Figure  19c  shows  that  the  average  improvement  factor  for  a  2-port  doppler  processor 
(dotted  curves)  and  the  improvement  factor  for  a  2-pulse  MTI  (dashed  curves)  are  identical 
unless  the  clutter  spectrum  is  very  wide.  The  reason  for  this  is  that  the  optimum  weights 
for  the  2-port  processor  are  unaffected  by  the  addition  of  white  noise  until  the  noise  level 
becomes  high  enough  to  overcome  the  effects  of  the  Gaussian  clutter  spectrum.  This  does 
not  happen  as  long  as  the  noise  level  is  less  than  the  clutter  level  ( Ew  <  0.5). 

Figures  20c,  21d,  and  22d  show  that  the  N- port  doppler  processor  for  N  >  2  gives  a 
better  average  improvement  factor  than  the  N-pulse  MTI  improvement  factor  even  when 
the  clutter  spectrum  is  narrow.  As  N  is  increased,  this  advantage  of  the  ,V-port  doppler 
processor  becomes  greater,  so  that  the  need  for  an  optimal  design  becomes  rac.e  important. 
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(b)  Weight  2,  Filter  1 


Fig.  12  —  Optimum  weights  for  a  2-port  processor.  □  represents 
the  real  part  of  the  complex  weights.  X  represents  the  imaginary 
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(d)  Weight  2,  Filter  2 


Pig  12  (Continued)  -  Optimum  weights  for  a  2-port  processor. 
,j  represents  the  real  part  of  the  complex  weights.  X  represents 
the  imaginary  part. 
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Fig.  13  —  Transfer  characteristics  of  the  2-port  doppler 
processor,  u  represents  both  filters  when  the  clutter  spectrum 
is  narrow  (o CT  <  0.15).  x  represents  filter  1  and  A  represents 
filter  2  when  ofT  >  0.15. 


Figure  18  shows  the  four  filter  shapes  of  a  4-port  doppler  processor  for  three  values 
of  clutter  spectral  widths  (acT  =  0.005,  0.1,  0.5).  If  the  clutter  is  narrow,  all  four  filters 
are  the  same  and  are  also  the  same  as  in  the  4-pulse  optimal  MTI.  If  the  clutter  is  wide, 
all  four  filters  approach  a  (sin  x)/x  shape  centered  at  their  respective  detection  regions. 

The  transition  of  these  filter  shapes  from  an  optimal  4-pulse  MTI  to  (sin  x)/x  filters 
centered  at  their  detection  regions  is  shown  in  Fig.  23.  If  white  noise  is  10'4  of  the  total 
interference,  these  filter  shapes  are  altered  appreciably,  as  shown  in  Fig.  24.  The  effect  of 
white  noise  is  much  more  severe  for  narrow  clutter.  For  wide  clutter  all  four  filters  again 
approach  (sin  x)/x. 

Comparing  Fig.  24(a)  with  Fig.  23(a),  we  see  that  filter  1  is  drastically  affected  by 
noise,  especially  for  a  narrow  clutter  spectrum.  With  the  addition  of  noise  the  optimum 
filter  is  no  longer  a  4-pulse  MTI;  it  is  seen  that  the  peak  response  is  moved  closer  to  the 
center  of  the  detection  region  (zero  doppler).  If  Fig.  24b,  c,  and  d  is  compared  with 
Fig.  23b,  c,  and  d,  similar  effects  can  be  seen  for  filters  2,  3,  and  4.  The  peak  response 
moves  closer  to  the  detection  region  and  the  sidelobe  level  increases  (and  the  main  lobe 
level  narrows)  for  a  narrow  clutter  spectrum.  These  effects  become  more  pronounced  as 
the  fractional  noise  energy  is  increased,  0.01  for  Fig.  25  and  0.5  for  Fig.  26. 

TTie  results  of  this  investigation  show  that  the  clutter-to-noise  ratio  is  an  important 
parameter  for  determining  the  optimal  iV-port  doppler  processor. 


CONCLUSIONS 

An  optimization  procedure  has  been  developed  for  the  transversal  filter  using  the 
maximum-likelihood-ratio  criterion.  When  the  interference  has  Gaussian  amplitude 
statistics,  the  likelihood  ratio  is  maximized  when  the  output  signal-to-interference  ratio  is 
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maximized.  Under  these  conditions,  the  optimum  weights  for  the  transversal  filter  can  be 
derived  from  the  covariance  of  the  signal  and  the  interference. 

If  Mj  is  the  covariance  matrix  of  the  interference  and  Ms  is  that  of  the  signal,  the 
maximum  improvement  in  signal-to-interference  ratio  (the  improvement  factor)  is  given  by 
the  largest  eigenvalue  of  M's  where  is  defined  by  Eq.  (34).  The  set  of  optimum 
weights  is  given  by  the  eigenvector  Mg  that  is  associated  with  this  largest  eigenvalue.  With 
the  design  procedure  reduced  to  finding  the  eigenvalues  and  eigenvectors  of  a  matrix,  the 
filter  for  a  given  application  can  be  designed  using  numerical  techniques  to  compute  the 
eigenvectors  and  eigenvalues.  The  optimal  design  can  be  computed  with  this  procedure 
and  the  improvement  factor  can  be  obtained  with  no  additional  computing. 

The  covariance  of  a  function  of  time  is  the  Fourier  transform  of  the  power  spectrum 
of  the  function.  Therefore,  the  spectrums  of  the  signal  and  the  interference  must  be 
known  prior  to  deriving  the  optimum  weights.  This  optimization  procedure  has  been 
applied  to  the  detection  of  a  signal  with  an  unknown  doppler  shift.  Three  cases  were 
considered  in  this  research.  These  were  defined  in  terms  of  the  a  priori  knowledge  of  the 
doppler  shift,  as  follows:  (a)  the  doppler  shift  is  completely  unknown,  (b)  the  doppler 
shift  is  known  exactly,  and  (c)  the  doppler  shift  is  known  within  a  region. 

It  was  shown  that  if  the  doppler  shift  is  completely  unknown  this  optimization 
procedure  gives  the  same  results  as  the  procedure  of  Emerson  (3),  which  is  to  maximize 
the  rejection  of  interference.  However,  it  has  been  shown  in  this  research  that  his  pro¬ 
cedure  is  equivalent  to  one  using  the  stronger  criterion  of  maximizing  the  output  signal- 
to-interference  ratio. 

The  procedure  developed  in  this  research  is  equivalent  to  that  of  Applebaum  (4)  and 
Brennan  (7)  when  the  doppler  shift  is  known  exactly.  For  this  condition  a  filter  is 
optimized  to  detect  a  signal  at  a  given  doppler  shift. 

When  the  doppler  shift  is  unknown,  it  must  be  assumed  to  have  equal  probability  of 
occurring  at  any  frequency  within  the  analyzing  bandwidth  of  the  filter.  However,  this 
bandwidth  can  be  covered  by  several  contiguous  filters  instead  of  a  single  filter.  This  in¬ 
curs  no  additional  data  storage  when  several  independent  sets  of  weights  are  applied  to  the 
data.  In  this  way,  each  filter  can  be  designed  to  detect  a  signal  within  a  portion  of  the 
total  bandwidth  instead  of  the  entire  bandwidth.  The  doppler  shift  of  the  signal  can  be 
assumed  to  have  equal  probability  of  occurring  anywhere  within  a  region  instead  of  the 
entire  doppler  domain.  This  additional  information  about  the  signal  has  been  used  to  im¬ 
prove  the  output  signal-to-interference  ratio. 

The  results  of  this  research  were  applied  to  a  radar  incorporating  a  single-port  MT1 
processor.  In  this,  a  single-filter  output  is  used  to  detect  returns  from  moving  targets 
which  may  have  any  doppler  shift  and  to  reject  returns  from  fixed  objects.  These  latter 
returns  are  referred  to  as  clutter.  A  zero-mean  Gaussian  clutter  spectrum  was  used  to 
compare  the  improvement  factor  of  the  optimal  MTI  with  that  of  the  conventional  MTI. 
The  results  are  shown  in  Fig.  8,  which  shows  that  the  increased  improvement  factor  i3 
affected  only  slightly  by  the  width  of  the  clutter  spectrum  or  standard  deviation  of  the 
Gaussian  spectrum.  Typical  values  for  this  increase  are  0  dB  for  a  2-pulse  MTI,  2  dB  for 
a  3-pulse  MTI,  3  dB  for  a  4-pulse  MTI,  and  5  dB  for  a  5-pulse  MTI.  The  advantage  of 
using  optimum  weights  is  modest  when  the  returns  from  only  a  small  number  of  pulses 
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Fig.  14  —  Optimum  weight*  fo"  filter  1  of  *  4-port  doppler 
processor.  □  represent*  the  re*l  p«rt  of  the  complex  weight*.  X 
represent*  the  imaginary  part. 
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(d) 

Fig.  14  (Continued)  —  Optimum  weights  for  filter  1  of  a  4-port 
doppler  processor.  D  represents  the  real  part  of  the  complex 
weights.  X  represents  the  imaginary  part. 
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Fig.  15  —  Optimum  weights  for  filter  2  of  a  4-port  doppler 
processor.  □  represents  the  real  part  of  the  complex  weights.  X 
represents  the  imaginary  part. 
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(d) 

Fitf.  15  (Continued)  —  Optimutn  weights  for  filter  2  of  a  4-port 
doppler  processor.  □  represents  the  real  part  of  the  complex 
weights.  X  represents  the  imaginary  part. 
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(b) 

Fig.  17  —  Optimum  weights  for  filler  '1  of  a  4 -port  dopplcr 
processor.  □  represents  the  real  part  of  the  complex  weights.  X 
represents  the  imaginary  part. 
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Fig.  17  (Continued)  —  Optimum  weights  for  filter  4  of  a  4-port 
floppier  processor.  1 1  represents  the  real  part  of  the  complex 
weights.  .V  represents  the  imaginary  part. 
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(c)  Average  (dotted  curve)  and  optimal  3-pulse 
MTI  (dashed  curve) 
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are  processed  at  one  time,  but  it  quickly  becomes  significant  as  the  number  of  pulses  for 
which  returns  are  processed  is  increased. 

The  Uchniques  developed  in  research  were  also  applied  to  a  radar  tf-port  doppler 
processor  in  which  N  independent  sets  of  filter  weights  were  used  to  design  N  filters.  For 
this,  the  doppler  domain  was  divided  into  N  regions.  Each  filter  was  optimized  for  a  tar¬ 
get  whose  doppler  shift  has  equal  probability  of  occurring  anywhere  within  one  of  the 
regions.  These  results  are  shown  in  Fig.  11.  Comparing  the  average  improvement  factor 
of  the  N  filters  of  an  N- port  processor  with  the  improvement  factor  of  the  optimal  MTI 
shows  that  the  optimal  MTI  approaches  a  0-dB  improvement  factor  when  the  clutter 
spectrum  is  wide.  The  IV-port  processor  provides  an  improvement  in  signal -to-clutter  ratio 
even  under  this  condition.  The  value  of  this  improvement  factor  depends  on  the  number 
of  returns  processed.  Typical  values  are  2.5  db  for  a  2-port  processor,  4  dB  for  a  3-port 
processor,  5  dB  for  a  4-port  processor,  and  6  dB  for  a  5-port  processor. 

The  optimum  weights,  or  filter  shapes,  vary  as  the  width  of  the  clutter  spectrum  is 
changed.  This  is  shown  in  Fig.  18  for  a  4-port  processor.  When  the  width  of  the  clutter 
spectrum  is  very  small,  all  four  filters  are  the  same.  Furthermore,  they  are  the  same  as 
the  optimal  MTI  with  a  single  output.  Therefore,  when  the  width  of  the  clutter  spectrum 
is  very  small,  only  one  output  is  needed,  and  that  output  is  the  optimal  MTI.  Unexpect¬ 
edly,  the  peak  response  of  the  filters  was  sometimes  found  to  occur  outside  the  detection 
region  for  that  particular  filter.  This  implies  that  for  clutter  with  a  very  narrow  spectrum, 
maximization  of  the  improvement  factor  results  mainly  from  minimizing  the  output 
dutter. 

For  clutter  with  a  large  spectral  width,  the  characteristic  of  each  of  the  four  filters 
identified  by  triangles  on  the  curves  is  very  close  to  a  (sin  x)/x  function  in  shape,  and  the 
peak  response  is  centered  in  the  detection  region  of  each  filter.  This  implies  that  for 
dutter  of  large  spectral  width,  the  optimal  processor  is  a  discrete  Fourier  transform  with 
no  weighting  of  the  input  data. 

In  the  transition  region  between  small  and  large  spectral  widths  of  the  clutter,  the 
peaks  of  the  filter  responses  are  closer  to  the  center  of  the  detection  regions.  In  this 
region,  the  sidelobes  of  the  filter  response  are  found  to  be  very  high,  even  higher  than 
those  of  the  (sin  x)/x  filter  shape.  The  shapes  of  the  filters  are  controlled  by  weighting 
the  input  data. 

Hie  optimum  filter  weights  are  easily  computed  when  the  clutter  has  either  a  very 
narrow  or  a  very  wide  spectrum.  In  the  transition  region,  it  is  necessary  to  go  through 
the  optimization  procedure  developed  during  this  research.  The  design  parameters  of  a 
radar  system  must  be  chosen  with  many  factors  in  mind.  In  general,  they  cannot  be 
selected  in  such  a  way  that  the  ratio  of  the  clutter  spectrum  and  the  PRF  is  either  very 
small  or  very  large.  For  most  applications,  the  doppler  processor  must  be  designed  to 
operate  within  this  transition  region.  Furthermore,  as  the  number  of  returns  to  be 
processed  is  increased,  this  transition  region  becomes  wider.  That  is,  the  clutter  spectrum 
must  be  much  narrower  before  all  the  filters  approach  the  optimal  MTI,  or  the  spectrum 
must  be  much  wider  before  the  IV-port  processor  approaches  the  discrete  Fourier  trans¬ 
form.  Therefore,  for  a  radar  system  processing  a  large  number  of  returns  this  optimiza¬ 
tion  procedure  is  necessary  to  ensure  the  best  performance  in  detecting  moving  targets. 
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(d)  Average  (dotted  curve)  and  optimal  4-puise  MTl  (daahed  curve) 

Fi«.  21  (Continued)  —  The  effect  of  white  noise  on  the  improvement  factor  of 
a  4-port  doppler  processor.  Ew  is  the  fraction  of  the  total  interference  energy 
that  ia  white  noise. 
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(a)  Filter  on* 


(b)  Filter  two,  fir* 

Fig.  22  —  The  effect  of  white  noise  on  the  improvement  factor  of  a  5-port  dop- 
pier  processor.  Ey/  is  the  fraction  of  the  total  interference  energy  that  is  white 
noiae. 


42 


(d)  Average  (dotted  curve*)  and  optimal  5-pulae  MTI  (dashed  curvet) 


Pi*.  22  (Continued)  —  The  effect  of  white  noiae  on  the  improvement  factor  of 
a  5-port  doppler  processor.  Ew  it  the  fraction  of  the  total  interference  energy 
that  is  white  noiae. 
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F>*.  23  —  The  effect  of  the  width  of  the  clutter  spectrum  (ocT)  on  the  niter  shapes 
of  a  4-port  doppler  processor.  No  white  noiae. 


0.(1 


^fj||i 'lUirpii^iy I1  n^i?FT^T]m^q!m^ripii p ir^iipn^npiii^^  O' 

DOPPLER  EREQUENCY/RRDfiR  PRF 


(c)  Filter  No.  3,  centered  at  1/2  PRF 


FV  23  (Continued)  —  The  effect  of  the  width  of  the  clutter  spectrum  (oc7^ 
the  filter  shape*  of  a  4-port  doppler  processor.  No  white  noise. 
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Fig.  24  (Continued)  —  The  effect  of  the  width  of  the  clutter  spectrum  (ocT) 
on  the  filter  shape*  of  a  4-port  doppler  processor.  The  white  noise  input  energy 
is  10-4  of  the  total  input  interference  energy. 
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Fig.  25  —  The  effect  of  "the  width  of  the  clutter  spectrum  ( acT )  on  the  filter  shapes 
of  a  4-port  doppler  processor.  The  white  noise  input  energy  is  10*2  of  the  total  in¬ 
put  interference  energy. 
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Fig.  25  (Continued)  -  The  effect  of  the  width  of  the  clutter  spectrum  ( o(T ) 
on  the  Filter  shapes  of  a  4-port  doppier  processor.  The  white  noise  input  energy 
is  10— *  of  the  total  input  interference  energy. 
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Fia  26  -  The  effect  of  the  width  of  the  clutter  spectrum  (acT)  on  the  filler  shapes 
oU  4  port?oppler  proces^r.  The  whit,  noise  input  energy  is  0.5  of  the  total  mput 
interference  energy. 
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Fig.  26  (Continued)  —  The  effect  of  the  width  of  the  clutter  spectrum  {acT) 
on  the  filter  shapes  of  a  4-port  doppler  processor.  The  white  noise  input  energy 
is  0.5  of  the  total  input  interference  energy. 
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It  has  also  been  shown  by  this  research  that  neither  the  conventional  nor  the  optimal 
MTI  provides  any  improvement  factor  when  the  interference  corresponds  to  white  noise. 
Therefore,  when  the  total  interference  consists  of  both  clutter  and  white  noise,  the  im¬ 
provement  factor  is  reduced  significantly  (Figs.  9  and  10).  The  A/ -port  processor  provides 
a  larger  improvement  factor  against  clutter  plus  white  noise  (Figs.  19-22).  The  shapes  for 
the  optimal  filters  of  a  4-port  processor  are  shown  in  Fig.  23-26  for  noise  levels  constitut¬ 
ing  fractions  (0,  0.0001,  0.01,  and  0.5)  of  the  total  interference.  As  the  noise  level  is 
increased,  the  optimal  filter  shapes  never  approach  the  optimal  MTI  for  narrow  clutter, 
even  when  the  ratio  of  the  standard  deviation  of  the  clutter  to  the  PRF  of  the  radar  is  as 
small  as  0.005.  However,  the  optimal  shape  approaches  that  of  the  (sin  x)/x  function 
when  the  spectral  width  of  the  clutter  is  smaller  than  necessary  for  interference  consisting 
only  of  clutter.  Therefore,  the  width  of  the  transition  region  seems  to  be  unaffected  by 
the  addition  of  white  noise.  The  transition  region  as  a  whole,  though,  is  shifted  down¬ 
ward. 

Since  the  shapes  of  the  optimal  filters  are  altered  drastically  by  die  addition  of  white 
noise,  the  clutter-to-noise  ratio  is  an  important  parameter  for  determining  the  optimal  N- 
port  doppler  processor.  Therefore,  the  optimization  procedure  developed  during  this  re¬ 
search  will  ensure  the  best  performance  in  the  most  important  and  general  case,  in  which 
a  radar  system  must  contend  with  both  clutter  and  white  noise. 
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Appendix  A 

THEORY  OF  TRANSVERSAL  FILTERS 


In  many  fields  of  science  and  technology,  signals  are  detected  or  analyzed  on  the 
basis  of  their  spectral  content.  In  such  fields  as  communications,  radar,  and  sonar,  the 
transmitted  signal  is  designed  to  have  special  spectral  characteristics.  These  characteristics 
are  used  at  the  receiver  to  detect  and  identify  the  signal.  In  such  fields  as  medical 
technology,  radiometry,  and  passive  sonar,  the  receiver  designer  does  not  have  control  of 
the  signal  characteristics;  however,  particular  spectral  characteristics  can  be  related  to 
particular  physical  occurrences.  The  fundamental  theory  used  in  these  and  many  other 
branches  of  science  and  technology  is  that  associated  with  linear  filtering  and  spectrum 
analysis.  A  review  of  digital  filtering  theory  can  be  obtained  from  the  work  of  Gold  and 
Rader*.  The  most  important  difference  between  digital  and  analog  filtering  is  that  in  the 
former  the  input  signal  must  be  sampled.  Assuming  a  constant  sampling  rate,  the  sam¬ 
pling  interval  determines  the  highest  frequency  that  can  be  analyzed.  Frequency  com¬ 
ponents  higher  than  this  maximum  are  “folded”  back  into  the  analyzing  bandwidth  and 
cause  “aliasing.” 

The  canonical  form  for  a  digital  filter  is  defined  as  that  form  which  requires  the 
minimum  storage,  or  memory,  to  perform  a  particular  operation.  From  Gold  and  Rader 
(Al),  the  canonical  form  for  an  Nth-order  filter  is  shown  in  Fig.  Al.  The  output  y(f)  for 
this  filter  can  be  described  in  terms  of  the  input  x(t)  by  an  Nth-order  difference  equation, 

N  N 

y(nT)  *  £  apinT-iT)-  £  b^nT-iT),  (Al) 

i-0  i-l 


where  T  is  the  sampling  interval.  The  input  and  output  are  defined  only  at  the  particular 
times  t  »  nT.  The  form  of  Eq.  (Al)  shows  the  recursive  relationship  i.eeded  to  derive 
output  y(t)  at  time  t  *  nT.  The  nonrecursive  form  of  Eq.  (Al)  is  obtained  by  setting 
bj  *  0,  for  all  i.  The  canonical  form  for  a  nonrecursive  digital  filter  is  shown  in  Fig.  A2. 
All  digital  filters  can  be  classed  as  either  recursive  or  nonrecursive. 


Transversal  Filters 

A  transversal  filter  is  a  nonrecursive  filter  as  shown  in  Fig.  A2.  The  input/output 
relationship  for  this  filter  is 


s 

y(nT )  -  22  ai*(nT  ~  «T).  (A2) 

i-0 


*B.  Gold  and  C.  M.  Rader,  Digital  Proceaai-.g  of  Signala,  McGraw-Hill,  New  York,  1969. 
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Pig.  A2  —  Canonical  form  for  an  nth-order  non  recursive  digital  filter 


This  equation  can  also  be  expressed  as  either  of  two  matrix  operations, 

y(nT)  -  aTx 

or 

y(nT)  *  xTa, 

where 


■x(nT)  ' 

*1  ' 

°2 

x(nT  -  T) 

*2 

a  * 

• 

and  x  ■ 

• 

m 

• 

?1 S*\ 

x(nT-NT) 

1 

Subscript  T  represents  the  transpose  of  the  matrix. 
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The  output  power,  as  shown  by  Applebaum* 

P  -  ly(nT)|2  -  yy^.. 

The  asterisk  represents  the  complex  conjugate.  If  Eqs.  (A3)  and  (A4)  are  used, 

?“«TxxJ.a*.  (A5) 

Expected  Output  Power— The  expected  output  power  is  found  by  taking  the  ex¬ 
pected  value  of  Eq.  (A5): 

P  ■  £  jaTxx5.a*j  *  aTxx^a*. 

If 

Af  -  xx*, 

then 


P-ayAfa*,  (A6) 

where  Af,  the  covariance  matrix  of  the  input,  can  be  derived  from  the  Fourier  transform 
of  the  normalized  input  power  density  spectrum.  The  normalization  is  such  that 

mm 

J  |x(f)|2dt  -  1 

or,  using  Parse val’s  theorem, 

■» 

J  Pc(/w- 1. 

where  pc {f)  -  |Jf(f)|2  and  X{f)  is  the  Fourier  transform  of  x(t). 

Power  Transfer  Function— The  power  transfer  function  is  usually  found  by  taking 
the  ratio  of  the  output  power  to  the  input  power  as  a  function  of  the  frequency  of  the 
input.  With  the  normalization  described  above,  the  power  of  a  single-frequency  input  is 
unity,  and  it  is  only  necessary  to  compute  the  power  at  the  output  as  a  function  of 
frequency.  Also,  since  the  filter  characteristic  repeats  at  multiples  of  the  sampling  fre¬ 
quency  (1  IT),  it  is  only  necessary  to  compute  the  output  power  for  frequencies  between 
zero  and  1  IT.  Consider  an  input: 

x(f)  -  cos  +  <p).  (A7) 


I 

i 

i 

i 


i 


.i 


•S.  P.  Applebaum,  “ Adaptive  Arrays,”  Syracuse  University  Research  Corp.  SPL-769,  June  1964. 
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If  N  is  the  number  of  delay  lines  (or  memory  elements)  in  the  transversal  filter  and  T  is 
the  time  interval  between  samples,  NT  is  the  total  processing  time  of  the  filter.  N  +  1 
samples  are  taken  during  this  time.  When  sampled  at  t  •  iT,  Eq.  (A7)  becomes 

N 

*(0*  2]  ®(f  “  'T>  co*  <2jrf0f  +  0). 

1*0 

The  covariance  of  this  input  is  found  by 

*(r)-  j  x(t)x(t-r)dt  (A8) 


'  N 

'  N 

T]  6(t  -kT)  cos  (2vfQt  +0) 

y  5(f  -  t  -  iT)  cos  (2wf0(t  -  r)  +  0) 

.Ar-0 

.  i-0 

The  integral  is  zero  except  when  t  -  kT  -  0  and  t-r-iT»0orr»(fc-  i)T.  Therefore 
the  (i,  k)  term  of  the  covariance  matrix  is 

mi  k  m  cos  (2 xf0kT  +  0)  cos  (2*f0iT  ♦  0),  (A9) 

where  i  -  0,  1,  ...,  N  and  k  *  0,  1,  ....  N.  To  normalize  M,  it  is  necessary  for 

j  x(t)2dt  ■*  j  x2(t)dt »  1 

When  t  ■  0,  Eq.  (A8)  becomes 

0(0) «  J  x2(t)dt. 

—mm 

But  r  ■  0  implies  k  ■  i  in  Eq.  (A9),  and  therefore  M  is  normalized  if 

0(0)  -  m..  -  cos2(2ff/0<T  +  0)  -  1. 

Since  phase  0  is  arbitrary  it  is  chosen  to  be  0  -  -2irf0iT.  Then  M  is  normalized  and  mi  k 
becomes 

mih  -  cos  2t rfQ(k  -  i)T.  (A10) 

In  summary,  the  power  transfer  function  of  a  transversal  filter  can  be  derived  by 
using  Eq.  (A6)  with  a  real  input  from  Eq.  (A7),  which  results  in  a  covariance  matrix  de¬ 
fined  by  Eq.  (A10). 

Complex  Filter  Weights—  For  most  applications  of  digital  filters,  it  is  desirable  for  the 
spectrum  being  filtered  to  be  at  “baseband.”;  that  is,  the  bandwidth  should  span  a  region 


/ 
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Fig.  A3  —  Transition  of  a  real  zpectrum  to  baseband, 
(a)  Real  lignal  to  be  analyzed,  (b)  Sinuaoidal  refer¬ 
ence  signal,  (c)  Spectrum  (a)  times  spectrum  (b). 
(d)  Aliasing  caused  by  using  a  reference  lignal  greater 
than  fg-8/2. 


f0 


(b) 


8/2  0  8/2 


Fig.  AS  —  Complex  translation  of  a  real  spectrum  to  baseband, 
(a)  Complex  exponential  reference  signal,  (b)  Spectrum  of  Fig. 
A3(a)  times  complex  reference  signal  (a). 


from  zero  to  some  maximum  frequency.  The  maximum  frequency  determines  the  mini¬ 
mum  sampling  frequency  to  meet  the  Nyquist  sampling  criterion.  If  the  signal  spectrum 
of  interest  is  not  at  baseband,  it  is  usually  multiplied  by  a  sinusoidal  signal  of  appropriate 
frequency  to  translate  the  spectrum  to  baseband.  This  is  shown  in  Fig.  A3a,  b,  and  c.  If 
die  spectrum  is  translated  more,  then  aliasing  occurs  as  negative  frequencies  interfere  with 
positive  frequencies,  as  shown  in  Fig.  A3d. 

This  aliasing  problem  could  be  overcome  by  multiplying  the  signal  of  interest  with  a 
complex  exponential  signal.  The  spectrum  of  this  complex  signal  is  not  an  even  function, 
and  positive  frequencies,  therefore,  can  be  distinquished  from  negative  frequencies  as 
shown  by  Fig.  A4. 

A  transversal  filter  can  be  designed  to  have  a  complex  impulse  response  simply  by 
using  complex  weights  In  this  way,  a  filter  can  be  designed  which  discriminates  be¬ 
tween  positive  and  negative  frequencies.  An  example  of  this  is  the  discrete  Fourier  trans¬ 
form,  which  can  be  considered  a  contiguous  bank  of  transversal  filters.  N  filters  are 
formed  when  N  samples  are  taken.  The  effective  weights  of  the  feth  filter  are 

afca  »  '  n  =  \,2 . N  (All) 
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pulTwe  have  ""  U*in8  **  Weight*  defmed  by  *»•  <A11)  “d  a  complex  signal  in- 


*(<)  -  ej2*fQt,  (A12) 

the  power  transfer  function  of  a  discrete  Fourier  transform  can  be  derived  using  Eq.  (A6) 
ine  autocorrelation  of  Eq.  (A12)  is  M 

ej2wf°T , 

which  results  in  a  covariance  matrix  defined  by 

o(«-*)r 

The  power  transfer  function  for  9  of  the  16  filters  generated  by  a  16-sample  discrete 
Founer  transform  is  shown  in  Fig.  A5.  The  filter  characteristics  shown  in  Fig.  A5  repeat 
with  multiples  of  the  sampling  frequency  f,.  Therefore  the  filter  response  to  negative 
frequencies  from  zero  to  -f,  is  simply  the  shapes  of  these  filters  translated  to  the  negative 
^quency  region.  In  this  way,  an  unambiguous  analyzing  frequency  bandwidth  can  be 
defined  either  from  -fJ2  to  ♦/,/ 2  or  from  0  to  f,.  These  filters  can  discriminate  between 
positive  and  negative  frequencies  over  the  region  from  -fJ2  to  +f  / 2. 

>tfan!Ve[S!1  mten  Can  h*  imp,emented  with  real  operations.  Consider  real 
input  signal  x(t)  which  has  a  spectrum  and  center  frequency  fc,  as  shown  in  Fig.  A4  The 
spectrum  is  translated  to  baseband,  so  that 

x’(t)  •  x(t)ej2’fo*. 

When  an  input  signal  is  sampled,  a  vector  can  be  defined. 


where 


xt-je{(i  -  l)T]eywo(,'1)r 

-  x[(i  -  1)T]  (cos  2 jrf0(i  -  1  )T  +j  sin  2 nf0(i  -  1)T] 
Mxf  +jxf. 
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Doppler  Filtering 

It  is  well  known  in  the  fields  of  optics  and  acoustics  as  well  as  electromagnetics,  that 
if  there  is  relative  motion  between  the  source  of  oscillation  and  the  observer,  an  apparent 
shift  in  frequency  will  result.  This  is  the  doppler  effect  that  is  put  to  use  in  many  fields, 
such  as  radio  astronomy,  coherent  communication,  radar,  and  sonar.  Each  field  has  a 
need  to  measure  accurately  the  doppler  caused  by  relative  motion.  The  doppler  shift  fd 
of  a  frequency  ft  is 
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Fig.  A5  (Continued)  —  Transfer  characteristics  of 
a  16-sample  discrete  Fourier  transform.  Of  the 
16  filters,  9  are  shown;  the  remaining  7  Alters  are 
the  mirror  images  of  Alters  2  through  8. 
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where  c  is  the  speed  of  propagation  and  v  is  the  relative  velocity  between  source  and 
observer. 

Complex  Demodulation— Consider  a  signal  at  the  source, 

u(t)  ■  m(t )  cos  2jt ftt,  (A14) 

where  m(t)  is  the  modulation  and  ft  is  the  carrier  frequency.  With  relative  velocity  v  be¬ 
tween  source  and  observer,  this  signal  appears  to  the  observer  to  be 

o(t )  -  m(t  +  td )  [cos  2 it(ft  +  fd  )(t  +  td )] 
where  td  is  the  propagation  time. 

For  many  applications  of  this  waveform,  such  as  in  radar,  sonar,  and  communica¬ 
tions,  the  modulation  bandwidth  is  small  as  compared  with  the  carrier  frequency.  In  this 
case,  the  modulation  spectrum  is  altered  very  little  as  compared  with  the  shift  in  fre¬ 
quency  of  the  carrier,  so  that  the  receiver  spectrum  can  be  considered  to  be  the  same 
shape  as  that  of  the  transmitted  spectrum  but  translated  by  the  carrier  doppler  shift.  The 
received  spectrum  is  approximately 

v'(t)  •  m(f  +  td)co»  *fd)t*<t>] 


where  <t>  -  2*(/,  +  fd  )td . 

The  attenuation  due  to  propagation  loss  has  been  ignored  since  it  does  not  affect 
the  design  of  the  optimum  receiver  for  detecting  this  signal.  As  described  above,  this 
signal  can  be  translated  to  baseband  by  forming  the  product  x(t )  *  v'(t)e^2nfit  and  filter¬ 
ing  out  the  difference  frequency, 

x(f)  »  m(f  +  td)  [cos  (2ir(/(  +  fd)t  +  0)  cos  2wftt+j  cos  (2ir(ft  +  fd)t  +  <t>)  sin  2xf,t] . 
After  filtering,  we  have 

*(f)  =  ^  m(t  +  td)  [cos  (2xfdt  +  sin  (2ir fdt  +  <t>)] 

+  (A15) 

The  demodulated  complex  signal,  Eq.  (A15),  has  been  translated  to  baseband  and  is 
now  of  a  form  to  be  processed  by  a  digital  doppler  filter. 

Pulse  Doppler  Radar- The  doppler  shift  created  by  a  moving  target  can  be  measured 
with  a  pulse  doppler  radar.  Although  many  transmitted  waveforms  are  used  for  this  pur¬ 
pose  in  radar,  most  of  them  can  be  related  to  conventional  scanning  pulsed  radars.  These 
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radars  transmit  a  carrier  frequency  that  is  gated  on  and  off  by  a  sequence  of  pulses.  The 
envelope  of  the  basic  pulse  is 


where  Ep  is  the  energy  per  pulse  and  AT  is  the  pulse  width.  The  transmitted  waveform 
is  therefore 


i>(t)  =*  e>2nfcl  2]  P(t-nT ), 
no¬ 
where  f  is  the  carrier  frequency  and  T  is  the  interpulse  period.  The  composite  returns 
from  these  transmitted  pulses  are  made  up  of  moving  target  returns  and  usually  much 
larger  returns  from  fixed  objects  (clutter).  The  objective  is  to  design  a  processor  that  can 
pass  the  moving  targets  and  reject  the  clutter  by  processing  the  returns  from  a  fixed  num¬ 
ber  N  of  pulses. 

If  a  moving  target  is  at  azimuth  a  with  respect  to  the  axis  of  the  antenna  pattern 
G(d)  and  if  the  antenna  scans  at  a  rate  of  u>a  rad/s,  the  detailed  model  of  the  signal  re¬ 
turn  for  the  conventional  scanning  pulsed  radar  is 

•V 

x(t)*  aG2(cjft  £  p(t  -  nT  -  td). 

nm0 

The  carrier  frequency  has  been  removed  at  the  receiver.  The  term  a  *  Ae**  represents  the 
unknown  amplitude  and  phase  of  the  carrier  signal  return;  <j2irfdt)  is  the  doppler  modulation 
due  to  the  target  motion.  The  delay  ta  corresponds  to  the  distance  from  the  radar  to  the 
target.  Signal  parameters  a,  a,  fd,  and  td  are  not  known  a  priori. 

If  it  is  assumed  that  the  modulation  of  the  returned  signal  by  the  antenna  pattern 
G(0)  is  small  and  that  the  received  complex  amplitude  a  is  not  significant  to  the  receiver 
design,  then  the  signal  to  be  processed  by  the  doppler  filters  becomes 

N 

x(t)  -  ej2nfd(t'td)  £  p(t-nT-td),  (A16) 

n*  0 

which  is  of  the  same  form  as  Eq.  (A15). 

The  time  duration  AT  of  the  modulating  pulse  is  usually  very  short  as  compared 
with  the  period  of  the  doppler  frequency.  Therefore,  the  demodulated  signal  is  essen¬ 
tially  a  sampled  version  of  the  doppler-shifted  signal  where  the  sampling  rate  is  the  radar 
PRF.  Thus,  digital  filtering  technology  is  directly  applicable  to  radar  doppler  processing. 

For  analysis,  the  pulse  duration  can  be  considered  to  be  zero,  and  the  input  signal  to  the 
doppler  filter  becomes 
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.V 

x(t)  -  ei2,f>tit'td)  ]T  S(t-nT-td).  (A17) 

n-0 

With  Eq.  (A17)  as  an  input,  a  transversal  filter  can  be  designed  to  detect  moving  targets 
and  reject  returns  from  fixed  objects  (clutter). 


Appendix  B 

OPTIMUM  WEIGHTS  FOR  MTI  RADAR  PROCESSORS 


Table  B1 

Two  Cancelers  (Three  Pulses) 


oeT 

Improvement  Factor 
(dB) 

Optimum  Weights 

Real 

Imagining 

0.0050 

64,8868 

•0.5002 

•0.0000 

1.0000 

0.0000 

•0.5002 

0.0000 

0.0056 

6278868 

*0.5003 

-0.0000 

1.0000 

o.onoo 

-0.5003 

0.0000 

".0063 

6078876 

-0.5004 

-o.oooo 

1.0000 

0.0000 

-0.5004 

o.oooo 

0.0071 

56,8881 

•0.5005 

-0.0000 

1.0000 

0.0000 

-0.5005 

0.0000 

0 . 3y 79 

56.8888 

-0.5006 

-0.0000 

1.0000 

0.0000 

-0.5006 

o.oooo 

0 .0069 

54,8897 

-0.5008 

-0.0000 

1.0000 

0.0000 

•0.5008 

0.0000 

0 .0100 

52,8909 

-0.5010 

-o .oooo 

1.0000 

0 .0000 

•0.5010 

0.0000 

o  .  0112 

50.8924 

-0.5012 

-o.oooo 

1.0000 

0.0000 

-0.5012 

0.0000 

0.0126 

48,8943 

-0.5016 

-o.oooo 

1.0000 

o.oooo 

•0.5016 

0.0000 

0.0141 

46,8966 

-0.5020 

-0.0000 

1.0000 

0 .  oooo 

-0.5020 

0.0000 

0.0158 

44 , 8996 

-0.5025 

-0.0000 

1.0000 

o.oooo 

-0.5025 

0.0000 

0 .0177 

42,9034 

•0,5031 

-0.0000 
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Table  Bl— Continued 


ocT 

Improvement  Factor 
(dB) 

Optimum  Weight* 

Real 

Imagining 

1.0000 

0.0000 

-0.5031 

0.0000 

0 .0199 

40,9081 

-0.5039 

-o.oooo 

1.0000 

0.0000 

-0,5039 

0.0000 

0.0223 

38 , 9141 

-0.5049 

-o.oooo 

1.0000 

0.0000 

-0.5049 

o.oooo 

n.025i 

36,9218 

-0.5061 

0.0000 

1.0000 

0.0000 

-0.5061 

-0.0000 

0.0281 

34,9314 

-0.5077 

0.0000 

1.0000 

0.0000 

-0.5077 

0.0000 

0.0315 

32,9437 

-0,5097 

o.oooo 

1.0000 

O.OOOO 

-0,5097 

0.0000 

0 .  nj54 

30.9594 

-0.5121 

0.0000 

1.0000 

0.0000 

-0.5121 

o.oooo 

n .0397 

28,9794 

-0.5152 

-0.0000 

1.0000 

o.oooo 

-0.5152 

0 . 0000 

0 . 0446 

27,0051 

-0.5190 

0.0000 

1.0000 

O.OflOO 

-0.5190 

0.0000 

0.0500 

25  i  0381 

•0.5236 

0  .  onoo 

1.0000 

0.0000 

-0.5236 

o.oooo 

0.0561 

23,0807 

-0.5294 

-0.0000 

1.0000 

0.0000 

-0.5294 

o.oooo 

0 . 0  0  ?  9 

21.1360 

-0.5365 

-o.oooo 

1.0000 

0.0000 

-0.5365 

o.oooo 

0.O706 

19,2080 

*0.5451 

0.0000 

1.0000 

0.0000 

-0.5451 

0.0000 

0 . 0  /92 

17,3022 

-0.5554 

0 . onoo 

1.0000 

0.0000 

-0.5554 

-0.0000 

o  .00H9 

15,4258 

-0.5677 

0.0000 

1.0000 

0.0000 

-0.5677 

o  .onoo 
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Table  Bl— Continued 


oeT 

Improvement  Factor 
(dB) 

Optimum  Weights 

Real 

Imagining 

n.  3998 

13,5883 

-0.5819 

-0.0000 

1.0000 

o.onoo 

•0.5919 

o.onoo 

0.1119 

11.8021 

•0.5980 

-0.0000 

1.0000 

o.onoo 

•0.5980 

o.onoo 

0.1256 

10.0824 

-0.6157 

-0 .0000 

1.0000 

o.onoo 

-0.6157 

o.onoo 

3.1409 

8,4477 

-0.634? 

-o.onoo 

1.0000 

0 .0000 

-0.6342 

0.0000 

0.1531 

6,9188 

•0.6525 

-0.0000 

1.0000 

0.0000 

-0.6525 

o.onoo 

0.1774 

5,5175 

-0.6694 

-o.oooo 

1.0000 

0.0000 

-0.6694 

0.0000 

0.1991 

4,2642 

-0.6836 

-0.0000 

1.0000 

0.0000 

-0,6836 

o.onoo 

3.?23 3 

3.1756 

-0,6942 

-0.0000 

1.0000 

0.0000 

•  0.6942 

0.2536 

2.2616 

•0.7011 

1 

1.0000 

-0.7011 

0.2812 

1,5246 

-0.7048 

1,0000 

0.0000 

-0.7048 

0.3155 

0.9586 

-0,7064 

1.0000 

o.oooo 

•0,7064 

0.0000 

0,5513 

•0.7070 

1.0000 

0.0000 

•0.7070 

0.0000 

0.3S72 

0 , 2819 

-0.7071 

-o.onoo 

1.0000 

-0,7071 

0.0000 

n.4*56 

0,1236 

-0.7071 

-0.0000 

1.0000 

o.onoo 

-0,7071 

0 . 0OQ0 

0 . 5000 

0,0444 

-0.7071 

1.0000 

0.0000 

-0.7071 

0.0000 

O.  A.  ANDREWS,  JR. 


Table  B2 

Three  Cancelers  (Four  Pulses) 


oeT 

Improvement  Factor 
(dB) 

Optimum  Weights 

Real 

Imagining 

0.0050 

9570728 

-0.3337 

0.0000 

1.0000 

0.0000 

•1.0000 

-o.oooo 

0,3337 

0.0000 

0.0056 

92,1967 

-0.3337 

-0.0000 

1.0000 

0.0000 

-1.0000 

0.0000 

0.3337 

0.0000 

0.0063 

89,4701 

-0.3339 

0.0000 

1.0000 

0.0000 

-1.0000 

0.0000 

0.3339 

0.0000 

0.0071 

86,4063 

0.3340 

o.oooo 

-1.0000 

0.0000 

1.0000 

o.oooo 

-0.3340 

-0.0000 

0.0079 

83,4174 

-0.3342 

-o.oooo 

1.0000 

0.0000 

•1.0000 

0.0000 

0.3342 

0.0000 

0.0089 

80,4001 

0.3344 

-0.0000 

-1.0000 

0.0000 

1.0000 

0.0000 

-0.3344 

-0.0000 

0.0100 

77,4124 

0.3346 

-0  •  0000 

-1.0000 

0.0000 

1.0000 

0.0000 

-0,3346 

-o.oooo 

0.0H2 

74,4093 

•0 . 3350 

0.0000 

1.0000 

0.0000 

-1.0000 

-0.0000 

0.3350 

0.0000 

0.0126 

71,4135 

0.3354 

-0.0000 

-1.0000 

0.0000 

1.0000 

0.0000 

-0.3354 

-0.0000 

0.0141 

66,4156 

0.3359 

-0.0000 

-1.0000 

0.0000 

1.0000 

0.0000 

•0,3359 

-0.0000 

0.0150 

65,4193 

0.3366 

-0.0000 

•1.0000 

Q.onoo 

1.0000 

o.oooo 

-0.3366 

-0.0000 
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Table  B2— Continued 


oeT 

Improvement  Factor 
(dB) 

Optimum  Weights 

Real 

Imagining 

0.0177 

62,4244 

0.3375 

-0.0000 

*1.0000 

Q.onoo 

l.OOOu 

0.0000 

-0.3375 

-0.0000 

C ,  0199 

59^4311 

0.3385 

-0.0000 

•1.0000 

0.0000 

1.0000 

0.0000 

•0.3385 

-o.onoo 

0.0223 

56,4392 

•0.3399 

0.0000 

1.0000 

0.0000 

•1.0000 

-0.0000 

0.3399 

0.0000 

0.0251 

53,4498 

•0.3416 

0.0000 

1.0000 

0.0000 

•1.0000 

-0.0000 

0.3416 

0.0000 

0.0281 

50,4632 

•0.3437 

o.onoo 

1.0000 

0.0000 

•1,0000 

-o.onoo 

0.3437 

0.0000 

0.0315 

47,4804 

0.3463 

-0.0000 

*1.0000 

0.0000 

1.0000 

0.0000 

•  0.3463 

-o.onoo 

0.0354 

44,5024 

-0.3497 

o.onon 

1.0000 

0.0000 

•  1.0000 

-0.0000 

0,3497 

0.0000 

0.0397 

41,5309 

•  0.3538 

0.0000 

1.0000 

o.onoo 

*1.0000 

-o.onoo 

0.3538 

0.0000 

0.0446 

38  75678 

•0,3590 

0.0000 

1.0000 

o.onoo 

•1.0000 

-o.onoo 

0.3590 

0.0000 

0.0500 

3576160 

0.3655 

-o.onoo 

-1,0000 

0.0000 

1,0000 

0.0000 

•0.3655 

-o.onoo 

0.0561 

32,6793 

0.3736 

-0.0000 

-1.0000 

o.onoo 

1.0000 

o.onoo 

•0 . 3736 

-o.onoo 
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Table  B2— Continued 


oeT 

Improvement  Factor 
(dB) 

Optimum  Weights 

Real 

Imagining 

0.2233 

378424 

0.6032 

-o.oooo 

•1.0000 

o.oooo 

1,0000 

0.0000 

•0.6032 

•0.0000 

0.2506 

2,6936 

0,6112 

-o.oooo 

•1.0000 

0.0000 

1,0000 

0.0000 

•0.6112 

-o.oooo 

0.2612 

1,7920 

0.6155 

-0.0000 

-1,0000 

o.oooo 

1.0000 

o.oooo 

•0.6155 

-0.0000 

0.3155 

1,1155 

0.6173 

-0.0000 

•1.0000 

0.0000 

1.0000 

0.0000 

-0,6173 

-0 . onoo 

0.3540 

0,6367 

-0.6179 

0.0000 

1.0000 

0.0000 

-1.0000 

-O.OOOO 

0.6179 

0.0000 

0.3972 

0,3240 

0.6180 

-0 . 0000 

-1.0000 

0.0000 

1.0000 

o.oooo 

-0,6180 

-0.0000 

0.4456 

0,1417 

0,6180 

-0.0000 

•1,0000 

0.0000 

1.0000 

o.oooo 

-0.6180 

-o.oooo 

0.5000 

0,0508 

•0,6180 

0.0000 

1.0000 

0.0000 

•1.0000 

-0.0000 

0.6180 

0.0000 
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Table  B3— Continued 


oeT 

Improvement  Factor 
<dB) 

Optimum  Weights 

Real 

Imagining 

0.0251 

68,9416 

0.1750 

0.0000 

-0.6746 

0.0000 

1.0000 

0,0000 

•  0 , 6748 

0,0000 

0.1750 

0.0000 

0.0281 

64,9587 

0.1771 

-0.0000 

•0.6769 

0.0000 

1.0000 

0.0000 

•0.6769 

-0.0000 

0.1771 

0.0000 

0.0315 

60,9817 

0,1799 

-0.0000 

•0.6795 

0.0000 

1.0000 

0 • 0000 

•0.6795 

-0.0000 

0.1799 

C.0O00 

0,0354 

57,0111 

0.1833 

-0.0000 

•0.6827 

O.OOOfl 

- - 

1.0000 

0.0000 

•0.6827 

-0.0000 

0.1833 

0.0000 

0.0397 

53,0496 

0.1877 

-0.0000 

•0.6866 

O.onoo 

1.0000 

0.0000 

•  0.6868 

-O.OOOfl 

0.1877 

0.0000 

0.0446 

49,1001 

0.1932 

-0.0000 

•  0.6918 

0.0000 

1.0000 

0.0000 

•  0.6918 

-0.0000 

0.1932 

o.oooo 

0.0500 

45 ,1669 

0.2002 

-0.0000 

•0.6979 

0.0000 

1.0000 

0.0000 

•  0.6979 

-0.0000 

0,2002 

0.0000 

0.0561 

41,2562 

0.2090 

-0.0000 

•0.7054 

0.0000 

1.0000 

O.OOOfl 

•0.7054 

-O.OOOfl 

0,2090 

O.OOOfl 
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Tmble  B3— Continued 


oeT 

Improvement  Factor 
(dB) 

Optimum  Weights 

Real 

Imagining 

0.0629 

37,3763 

0,2200 

-o.ooon 

•0.7146 

0.0000 

1.0000 

0  .oooo 

•  0,7146 

-0.0000 

0,2200 

0.0000 

0.0706 

33,5394 

0.2338 

-0.0000 

•  0.7255 

Q.onoo 

1.0000 

0 . 0000 

•  0.7255 

-o.ooon 

0.2338 

0.0000 

0.0792 

29,7622 

0.2510 

-o.onoo 

•0.7304 

o.onoo 

1.0000 

0.0000 

•  0.7384 

-o.onoo 

0  •  2510 

0.0000 

0.0009 

26,0680 

0.2720 

-0.0000 

•0.7533 

0.0000 

1.0000 

0.0000 

•0.7533 

-0.0000 

0,2720 

o.onoo 

0.0998 

22,4001 

0.2971 

-0.0000 

•  0.7699 

o.ooon 

1.0000 

o.onoo 

•  0.7699 

-0.0000 

0.2971 

o.onoo 

0.1119 

19,0624 

0.3260 

-o.onoo 

•  0.7870 

o.onoo 

1.0000 

o.onoo 

•0.7870 

-o.onoo 

0.3260 

o.onoo 

0.1256 

15,8386 

0.3579 

-o.onoo 

•  0.8055 

o.onoo 

1,0000 

o.onoo 

•0.8055 

-o.onoo 

0.3579 

o.onoo 

0.1409 

12,8681 

0,3908 

-o.onoo 

•0.8223 

o.onoo 

1.0000 

O.ooon 

•0,8223 

-0.0000 

0.3908 

o.onoo 
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Table  B3— Continued 


oeT 

Improvement  Factor 
(dB) 

Optimum  Weights 

Real 

Imagining 

0.1581 

10,1984 

0.4219 

•o.onoo 

•0.8368 

0.0000 

1.0000 

0.0000 

•0.8368 

-o.onoo 

0.4219 

o.onoo 

0.1774 

7,8650 

0.4489 

-o.onoo 

•0.8481 

o.onoo 

1.0000 

o.onoo 

•0.8481 

•o.onon 

0.4489 

o.onoo 

0,1991 

5,8846 

0.4697 

-O.onon 

•0.8561 

O.onon 

1.0000 

O.onon 

•0,8561 

-o.onon 

0.4697 

o.onoo 

0.2233 

4,2538 

0.4841 

-o.onoo 

•0.8611 

o.onon 

1.0000 

o.onon 

•  0.8611 

-o.onon 

0.4841 

o.onon 

0.2506 

2,9522 

0.4928 

-o.onoo 

•0.8639 

o.onon 

1.0000 

o.onon 

•0.8639 

-O.onon 

0.4928 

O.onon 

0.2812 

1,9484 

0.4973 

-o.onon 

•0.8652 

o.onoo 

1.0000 

o.onoo 

•0.8652 

-o.onon 

0.4973 

o.onoo 

0.3155 

1,2056 

0.4992 

-O.onon 

•0.8658 

o.onon 

1.0000 

O.onon 

■0.6636 

-O.onon 

0.4992 

O.onon 

0.3540 

0,6853 

0.4998 

-O.onon 

•0.8660 

O.onoo 

1.0000 

O.onon 

•0.8660 

-O.onon 

0,4998 

O.onon 

0.3972 

0,3478 

0.5000 

-o.onoo 

•0.6660 

O.onon 

1.0000 

O.onon 

•0.8660 

-o.onoo 

0.5000 

o.onon 
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Table  B3— Continued 


oeT 

Improvement  Factor 
(dB) 

Optimum  Weights 

Real 

Imagining 

0.4456 

0,1519 

0,5000 

>0,0000 

■0,8660 

0.0000 

1.0000 

0.0000 

•0,8660 

-o.oooo 

0.5000 

0.0000 

0,5000 

O', 0544 

0.5000 

-o.onon 

■0,8660 

0.0000 

1.0000 

0.0000 

•0.8660 

>0.0000 

0.5000 

0.0000 

i 
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Appendix  C 

FILTER  CHARACTERISTICS  OF  N'-PORT  PROCESSORS 


BEST  AVAILABLE  COPT 


Fig.  Cl  —  Filter  transfer  characteristics  of  a  3-port  doppler 
processor.  J  represents  the  filter  shape  when  the  width  of  the 
clutter  spectrum  is  narrow  (0CT  *  0.005),  A  represents  the 
filter  shape  for  ocT-  0.1,  and  X  represents  the  filter  shape  for 
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